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Abstract 

Cortical neurons are subject to sustained and irregular synaptic activity which causes important 
fluctuations of the membrane potential (Vm)- We review here different methods to characterize this 
activity and its impact on spike generation. The simplified, fluctuating point-conductance model of 
synaptic activity provides the starting point of a variety of methods for the analysis of intracellular 

recordings. In this model, the synaptic excitatory and inhibitory conductances are described by 
Gaussian-distributed stochastic variables, or "colored conductance noise". The matching of experimen- 
tally recorded distributions to an invertible theoretical expression derived from the model allows 
the extraction of parameters characterizing the synaptic conductance distributions. This analysis can 
be complemented by the matching of experimental power spectral densities (PSDs) to a theoretical 
template, even though the unexpected scaling properties of experimental PSDs limit the precision of 
this latter approach. Building on this stochastic characterization of synaptic activity, we also propose 
methods to qualitatively and quantitatively evaluate spike-triggered averages of synaptic time-courses 
preceding spikes. This analysis points to an essential role for synaptic conductance variance in deter- 
mining spike times. The presented methods are evaluated using controlled conductance injection in 
cortical neurons in vitro with the dynamic-clamp technique. We review their applications to the analysis 
of in vivo intracellular recordings in cat association cortex, which suggest a predominant role for inhi- 
bition in determining both sub- and supra-threshold dynamics of cortical neurons embedded in active 
networks. 
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Cerebral cortical networks can generate states of intense and irregular activity, which are characterized by 
low-amplitude "desynchronized" fast activity in the electroencephalogram (EEG), a defining feature of the 
awake state. Intracellular measurements in awake animals (Woody and Gruen, 1978; Matsumura et al., 
1988; Baranyi et al., 1993; Steriade et al, 2001; Timofeev et al., 2001; Rudolph et al, 2007) have shown 
that cortical neurons are depolarized (about -60 mV), have a low input resistance, their membrane potential 
(Vm) fluctuates, and they fire irregularly and sustainedly. During slow-wave sleep, or under several types 
of anesthetics (such as urethane or ketamine-xylazine), the \m displays "up-" (depolarized) and "down-" 
(hyperpolarized) states, which are paralleled with EEG slow waves (Metherate and Ashe, 1993; Steriade et 
al., 1993; Steriade et al., 2001; Timofeev et al., 2001). During the up-state, the EEG is desynchronized and 
the V„j of cortical neurons is depolarized and highly fluctuating, similar to the sustained activity found in 
awake animals (Destexhe et al., 2007). Up- and down-states have also been found in ferret cortical slices 
using high-potassium and low-calcium extracellular media (Sanchez- Vives and McCormick, 2000) and in 
rat entorhinal slices as a function of the metabolic state (Cunningham et al., 2006): these experiments 
indicate that intracortical circuits are able to generate such states, presumably through recurrent excitation 
and inhibition. TTX block of action potentials in vivo (Pare et al., 1998) and CNQX block of excitatory 
synapses in vitro (Cunningham et al., 2006) abolish the depolarized, fluctuating states, which confirms their 
synaptic origin. 

The recurrent activity of cortical networks has been investigated using computational models at different 
levels and including various degrees of biological detail. Large networks of formal (or simple spiking) 
neurons allow the analytic derivation (or numerical confirmation) of conditions for different classes of 
network activity, such as oscillations or deterministic chaos (Van Vreeswijk and Sompolinsky, 1996; Roxin 
et al., 2005; Barak and Tsodyks, 2007). Network models incorporating a realistic diversity of cell types and 
details of cortical connectivity allow quantitative predictions to be obtained through numerical simulations 
(Compte et al., 2003; Hill and Tononi, 2005). At the single neuron level, detailed models can include 
the cell's morphology, a variety of intrinsic ion channels and distributed synaptic inputs: studies of such 
models assess the impact of massive input from the cortical network on dendritic processing, spike train 
statistics, neuronal responsiveness and integrative properties (Bemander et al., 1991; Destexhe and Pare, 
1999; Rudolph and Destexhe 2003b,c; Destexhe et al., 2003b). However, the parameterization of all those 
models requires a large amount of information that cannot be obtained from any single experiment, which 
complicates the comparison of simulation results with biological data. 

A complementary approach consists in developing and studying models that are relatively simple, but con- 
tain variables and parameters that can be related to quantities directly measured in experiments. Various 
mathematical studies of the firing dynamics of neurons with conductance-based synaptic inputs have been 
proposed (see for example Burkitt et al., 2003; Moreno-Bote and Parga, 2005; MuUer et al., 2007). Other 
authors investigated the integration of multiple current- or conductance-based inputs driven by Poisson 
spike trains, in terms of V^ fluctuations and output firing rate of simple point model neurons (Kuhn et 
al., 2004). In the approach proposed by Destexhe and colleagues (Destexhe et al., 2001), the effective im- 
pact at the soma of thousands of single synapses, respectively excitatory and inhibitory, activated by such 
Poisson spike trains, is represented by only two stochastic conductance variables, geit) and gi{t), mod- 
eled as Ornstein-Uhlenbeck (Brownian-motion-like) processes. This fluctuating point-conductance model 
describes the evolution of the subthreshold V^ of a single-compartment neuron by a passive membrane 
equation with the two additional stochastic conductance variables ge{t) and gi{t). We present here how this 
model can be used in close combination with electrophysiological experiments to investigate the properties 
and the impact of cortical recurrent activity at the single neuron level. 
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More specifically, we critically review, on the basis of new as well as already published data, applications 
of the point-conductance model of synaptic activity to the analysis of intracellular recordings in cortical 
neurons: in each case, we first briefly describe the method of analysis, then we show how it can be validated 
experimentally by using controlled "conductance injection" in biological neurons with dynamic-clamp, and 
finally we present the application of the method to the analysis of real synaptic activity. We especially focus 
on recently developed approaches for determining which of the conductance configurations occurring in 
the fluctuating synaptic activity trigger spikes. 



2 Methods 



2.1 Computational methods 

Computational models were based on single-compartment neurons described by the following membrane 
equation: 

dV 

C— = -GL{V-EL)-g,{V-E,)-gi{V-E,) + I,,t, (1) 

where C denotes the membrane capacitance, I^xt a stimulation current, Gl the leak conductance and the 
leak reversal potential. ge{t) and gi{t) are stochastic excitatory and inhibitory conductances, with respective 
reversal potentials Eg and Ei. 

These effective synaptic conductances were described by the following Omstein-Uhlenbeck model (Des- 
texhe et al., 2001): 

- 'fc(')-«.o] + ,/p4,W (2) 



dt X, 



e 



where ^^o and are, respectively, the mean value and variance of the excitatory conductance, Xg is the 
excitatory time constant, and ^e(0 is a Gaussian white noise source with zero mean and unit standard 
deviation. The inhibitory conductance gi{t) is described by an equivalent equation (Eq.|3]) with parameters 
gio, oj, Xi and noise source ^i{t). 

In some simulations, the voltage-dependent conductances responsible for action potentials, ^a^^, and g^d 
(with respective reversals E^a and Ek), were included. They were described by Hodgkin-Huxley type 
models (with equations and parameters identical as described in Destexhe et al., 2001, and references 
therein), resulting in the following equation: 

dV 

C— = -GL{V-EL)-gNaiV-ENa)-gKd{V-EK)-ge{V-Ee)-gi{V-E,) + Iext. (4) 

dt 



Simulations were performed on LINUX workstations using the NEURON simulation environment (Hines 
and Carnevale, 1997). 
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In vitro experiments were performed on 380-400 /jm thick coronal or sagittal slices from the lateral portions 
of pentobarbital-anesthetized adult ferret (Marshall Europe, Lyon) and guinea-pig (CPA, Olivet, France) 
occipital cortex, as described previously (Rudolph et al., 2004; Pospischil et al., 2007). Slices were main- 
tained in an interface style recording chamber at 33-35 °C in slice solution containing (in mM) 124 NaCl, 
2.5 KCl, 1.2 MgS04, 1.25 NaHP04, 2 CaCb, 26 NaHCOs, and 10 dextrose and aerated with 95% O2- 
5% CO2 to a final pH of 7.4. In some experiments on ferret cortical slices, after approximately 1 hour, 
the solution was modified to contain 1 mM MgS04, 1 mM CaCl2 and 3.5 mM KCl (Sanchez- Vives and 
McCormick, 2000) in order to obtain up-states. Intracellular recordings following two hours of recovery 
were performed in all cortical layers on electrophysiologically identified regular spiking and intrinsically 
bursting cells. 

All research procedures concerning the experimental animals and their care adhered to the American Phys- 
iological Society's Guiding Principles in the Care and Use of Animals, to the European Council Directive 
86/609/EEC and to European Treaties series no. 123, and was also approved by the local ethics committee 
"Ile-de-France Sud" (certificate no. 05-003). 

We also review data from intracellular recordings in cat association cortex in vivo, which were described 
in detail elsewhere (Rudolph et al., 2005; Steriade et al., 2001; Rudolph et al., 2007). 

2.3 Electrophysiology 

Sharp electrodes for intracellular recordings were made on a Sutter Instruments P-87 micropipette puller 
from medium-walled glass (WPI, IBFIOO) and beveled on a Sutter Instruments beveler (BV-IOM). Mi- 
cropipettes were filled with 1 .2-2 M potassium acetate - 4 mM potassium chloride and had resistances of 
65-1 10 MQ. after beveling. An Axoclamp 2B amplifier (Axon Instruments) was used for Ym recording and 
current injection. A Digidata 1322A card (Axon Instruments) was used for data acquisition at 20 kHz. 

The dynamic-clamp technique (Robinson et al., 1993; Sharp et al., 1993) was used to inject computer- 
generated conductances in real neurons. Dynamic-clamp experiments were run as described previously 
(Rudolph et al., 2004, Pospischil et al., 2007) using the hybrid RT-NEURON environment (developed by 
G. Le Masson, INSERM 358, Universite Bordeaux 2), which is a modified version of NEURON (Hines 
and Camevale, 1997) running in real-time under the Windows 2000 operating system (Microsoft Corp.). In 
these experiments, the injected current loynCiamp was determined from the fluctuating conductances ge{t) 
and gi[t) modeled with Ornstein-Uhlenbeck processes (Eqs. |2]-[3]) as well as from the difference between 
the recorded membrane voltage V and the respective reversal potentials: 

hynClamp = -geiV - Eg) - gi(V - Ei) . (5) 

The contamination of measured membrane voltage by electrode artefacts was avoided either through the 
use of the discontinuous current-clamp mode (in which current injection and voltage recording alternate 
at frequencies of 2-3 kHz with our electrodes) or with Active Electrode Compensation, a novel, high- 
resolution digital on-line compensation technique we have recently developed (Brette et al., 2005, 2007; 
Rudolph et al., 2005). 
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The PC-based software ELPHY (developed by G. Sadoc, CNRS Gif-sur-Yvette, ANVAR and Biologic), 
Statview, Excel, custom-written C-code and Neuron-code were used for analyses. All values are given as 
average ± standard deviation. 

2.4.1 VmD analysis 

In dynamic-clamp experiments re-creating up-states with conductance injection, two approaches were used. 
In 3 cells, conductance estimates were computed for different, realistic values of leak conductance and cell 
capacitance, and were then tested against the real up-states in dynamic-clamp. In 2 cells, leak conductance 
and cell capacitance were estimated from the response to short current pulses, and those values were then 
used for conductance estimation: the estimated conductance parameters, when used for dynamic-clamp 
injection, proved to allow successful matching of the real up-state in terms of distributions. In all cases, 
conductance estimation was done during the recording using ELPHY. 

VmD analysis of the in vivo data reviewed here is described in Rudolph et al., 2005 and Rudolph et al., 
2007. 

2.4.2 Power Spectral Density (PSD) analysis 

Power spectra of activity were fit to an analytic template (see Results) using a simplex fitting algorithm 
(Press et al., 1986). Different initial conditions ("first guesses") were given to the fitting procedure to 
ensure that there was no convergence to local minima. Some fits were realized by fitting both the amplitude 
of excitatory and inhibitory components (Ag and A,; see Results), as well as the time constants (Xg and T/). 
In other cases, it was not possible to fit 4 parameters from the experimental PSD. In such cases (typically 
from in vivo data), the fit was performed with a single amplitude component (Ag = A;). In all cases, the 
effective membrane time constant (Xm) was fixed to the value estimated from the recordings. 

2.4.3 Spike-Triggered Average (STA) analysis 

For each conductance injection, spikes were detected using a threshold at -30 mV. Inter-spike-intervals 
(ISIs) were computed and for further analysis, a stable region in terms of ISI distribution was used (as 
assessed by a non-significant Spearman correlation test between ISI duration and time), including 673 ± 
493 spikes. When investigating the impact of frequency on the accuracy of the estimates, the same number 
of spikes (52) was used for all analyzed injections. Spikes following ISIs of at least 100 ms were then 
selected. 50 ms-long pieces of Ym, ge and gi preceding each selected spike were averaged to obtain the 
"measured STAs". Using the recorded conductance traces, we checked that excluding spikes following ISIs 
shorter than 100 ms did not affect the measured STAs in an important way (the difference in conductance 
variation before the spike was of -0.14 ± 0.35 nS for excitation, and of 0.5 ± 0.7 nS for inhibition). 
However, in order to compare with the STA extraction method based on the STA, it is important to 
exclude short ISIs from the analysis to minimize contamination of the STA by conductances related to 
preceding spikes. 

For conductance STAs extraction based on the Y^. STA, parameters were obtained in the following way: 
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responses to depolarizing current pulses were used to estimate the membrane capacitance from the time 
constant of exponential fits to the decay of the V^, and the leak conductance was obtained by dividing the 
average voltage during conductance injection, relative to rest, by the average injected current. This results in 
an effective leak conductance possibly comprising different types of voltage-dependent conductances. An 
indication of the validity of this approach is provided by the estimated conductance STAs: the important 
variation of this effective leak conductance on the time scale of the length of the analysis window (for 
example due to conductances underlying spike frequency adaptation) would result in alterations of the 
stationary phase, which we did not observe when excluding ISIs shorter than 100 ms. More details about 
the application of the STA analysis to in vivo data can be found in Rudolph et al. (2007). 

To quantify both the measured and the extracted conductance STAs, we fitted the conductance time courses 
using the exponential template 



8eit) 



8eO 



l+ke exp 



jt-tp) 

rSTA 



(6) 



for excitation, and an equivalent equation for inhibition. Here, tp stands for the time of the spike, ke quan- 
tifies the maximal increase/decrease of conductance prior to the spike (Agg = ggQ kg), with time constant 



T^J^, and g^Q^^ is the average baseline conductance (see Results). 



,STA 



3 Results 

The point-conductance model of recurrent cortical activity describes the evolution of the subthreshold 

of a point neuron based on two effective fluctuating conductances gg and gt (Destexhe et al., 2001; see 
Methods for equations). The basic assumption of this model is that the synaptic conductances can be de- 
scribed as Gaussian-distributed stochastic variables. It was shown that the well-known Omstein-Uhlenbeck 
model of Brownian noise (Uhlenbeck and Omstein, 1930) approximates very well the total synaptic con- 
ductances resulting from a large number of simulated conductance-based synaptic inputs (Destexhe et al., 
2001; Destexhe and Rudolph, 2004). The fitting of a Gaussian model to the total synaptic conductances 
seen at the soma, in simulations performed using realistic cortical neuron morphologies and distributed 
synaptic inputs, indicates the following correspondences between variables: the time constants (Xg, x,) are 
identical to the decay time constants of synaptic currents. The average conductance (geo, gm) is related to 
the overall (integrated) conductance, which depends on the release frequency of the corresponding Poisson 
inputs, the quantal conductance and the decay time of synaptic currents. The variance of the conductances 
(Og, oj) is related to the same parameters, as well as to the amount of correlation between inputs of the 
same type (Destexhe et al., 2001). Correlations between presynaptic spike trains seem necessary in order 
to account for the high amplitude of fluctuations observed in vivo (e.g., Destexhe and Pare, 1999; Leger 
et al., 2005). 

It has to be noted that the Gaussian model of synaptic conductances can only be considered as an approxi- 
mation. As shown theoretically using the shot noise formalism, in the presence of very strong correlations, 
or a low mean presynaptic rate, this approximation might fail (Richardson and Gerstner, 2005; Rudolph and 
Destexhe, 2006; see also Kuhn et al., 2003). In this case, effective synaptic conductances might be better 
described by a gamma distribution. However, the Gaussian model seems appropriate, even in the presence 
of correlations, as long as the ratio of conductance standard-deviation to conductance mean remains small 
(so that the presence of negative conductances is negligible). Experimental support for this assumption has 
been obtained using voltage-clamp conductance measurements (Destexhe et al., 2003a). 
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The theoretical and numerical analysis of the point-conductance model has led to various useful derivations: 
an invertible expression for the steady-state distribution of V„j fluctuations (Rudolph and Destexhe, 2003a, 
2005), an expression for the power spectral density of V,„ fluctuations (Destexhe and Rudolph, 2004), a 
geometrical analysis of the configuration of conductances directly preceding spikes, and a probabilistic 
method for calculating the most likely conductance time course preceding spikes given an average 
time course (Pospischil et al., 2007). We now proceed to examine in more detail the applications of these 
computational results to the analysis of intracellular recordings in cortical neurons. 



3.1 The VmD method for extracting synaptic conductance parameters 
3.1.1 Outline of the VmD method 



The model described by Eqs. [T]-[3]has been thoroughly studied theoretically and numerically. This model 
describes the subthreshold V,„ fluctuations of a neuron subject to fluctuating conductances ge and gi. Dif- 
ferent analytic approximations have been proposed to describe the steady-state distribution of these Vm 
fluctuations (Rudolph and Destexhe, 2003a, 2005; Richardson, 2004; Lindner and Longtin, 2006; for a 
comparative study, see Rudolph and Destexhe, 2006). One of these expressions is invertible (Rudolph and 
Destexhe, 2003a, 2005), which enables one to directly estimate the parameters (geo, geo, <3e, CJ,) from ex- 
perimentally calculated Y,n distributions. This constitutes the basis of the VmD method (Rudolph et al., 
2004), which we outline below. 

The essential idea behind the VmD method is to fit an analytic expression to the steady-state subthreshold 
Vm distribution obtained experimentally, and yield estimates of the parameters (mean, variance) of the 
underlying synaptic conductances. Among the different analytic expressions outlined above, we consider 
the following Gaussian approximation of the steady-state V,,, distribution: 



p(y) ~ exp ^ ^ 



2o2 



(7) 



where V is the average V„, and Oy its standard deviation. This expression provides an excellent approxima- 
tion of the Ym distributions obtained from models and experiments (Rudolph et al., 2004), because the Vm 
distributions obtained experimentally show little asymmetry (for up-states and activated states; for specific 
examples, see Figs. [1] and |2| and Rudolph et al., 2004, 2005, 2007). 

One main advantage of this Gaussian approximation is that it can be inverted, which leads to expressions 
of the synaptic noise parameters as a function of the V,„ measurements, V and Oy. By fixing the values of 
Xe and X/, which are related to the decay time of synaptic currents and can be estimated from voltage-clamp 



data and/or current-clamp by using power spectral analysis (see Section 3.2), we remain with four param- 
eters to estimate: the means {geO, gio) and standard deviations (o^, o,) of excitatory and inhibitory synaptic 
conductances. To extract these four conductance parameters from the membrane probability distribution, 
Eq. |7]is, however, insufficient because it is characterized by only two parameters (V, Oy). To solve this 
problem, one possibility is to consider two V„, distributions obtained at two different constant levels of in- 
jected current Igxti and Igj^ti- In this case, the Gaussian approximation (Eq.|7]) of the two distributions gives 
two mean Vm values, Vi and V2, and two standard deviation values, Oyi and Oyi- The resulting system of 
four equations relating Vm parameters with conductance parameters can now be solved for four unknowns: 
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S{e,i}0 







{E{i,e} - 


-V2f 




{Ee - Vi ) {E, -V2) + iE,-V2) (E, - Vi ) 


{E{e,i}-E{ 


i,e})iVl-V2f 



{Ie.n -4xr2) - V2) + [Ie.,2 - Gl{E{,^,} - Ei)] {V, - V^) 

(%0-%4)(n-^2) 



^'''^ [{Ee - Vi ) {E, - V2) + {Ee - %) {E, - Vi ) J - (Vi - V2r 

Here, X{g ,} are effective time constants given by (Rudolph and Destexhe, 2005): 



where x„ 



'^{e,/} Xm 

C/ {Gl + + .?!o) is the effective membrane time constant. 



(8) 



(9) 



(10) 



These relations enable us to estimate global characteristics of network activity, such as mean excitatory 
{geo) and inhibitory (^,o) synaptic conductances, as well as their respective variances (a^, a?), from the 
sole knowledge of the V,„ distributions obtained at two different levels of injected current. This VmD 
method was tested using computational models and dynamic-clamp experiments (Rudolph et al., 2004) 
and was also used to extract conductances from different experimental conditions in vivo (Rudolph et al., 
2005, 2007; Zou et al., 2005). 



3.1.2 Testing the VmD method with dynamic-clamp 

Taking advantage of the possibility, given by the dynamic-clamp technique (see Methods), to mimic in a 
finely controlled way the fluctuating conductances ge and gi in biological neurons, we performed in vitro 
tests of the VmD method (Rudolph et al., 2004; Piwkowska et al., 2004). In a first test (in 5 neurons), 
we computed V^ distributions selectively during periods of subthreshold activity collected within up-states 
recorded in ferret cortical slices, we subsequently extracted conductance parameters from Gaussian fits 
to these distributions, and finally we used the estimated parameters to inject fluctuating conductances in 
dynamic-clamp in the same cell, during down-states. Fig [Ip shows a typical example of a real up-state 
and, shortly after, an up-state re-created in dynamic-clamp. We confirmed that the V,,, distributions are 
very similar in the two cases (see Rudolph et al., 2004 for more details). This test shows that the V,,, 
distributions observed experimentally in vitro during recurrent cortical activity can be accounted for by 
the proposed point-conductance model. We also re-estimated known parameters of synaptic conductances 
(ggO, giQ, CJg, Oi) injected in dynamic-clamp from the resulting V,„ distributions: the match between actual 
and estimated values is shown in Fig. [Ip. This second test indicates that the passive approximation for the 
membrane behavior holds in the studied case. In these tests, we did not consider the issue of the estimation 
of Xe and x,- and assumed these values are known. 



3.1.3 Analysis of intracellular recordings of cortical neurons in vivo 



The VmD method was then applied to analyze intracellular recordings in anesthetized (Rudolph et al., 
2005), as well as naturally sleeping and awake cats (Rudolph et al., 2007). 



Piwkowska et ah, J. Neurosci. Meth. (in press, 2008) 



9 



In the first study, recordings were performed in cat association cortex under ketamine-xylazine anesthesia, 
during the slow oscillation typical of this anesthetic and resembling slow-wave-sleep, as well as during 
prolonged periods of activity, triggered by brain stem (PPT) stimulation, and with activity similar to that 
of the aroused brain. The VmD method was used to extract synaptic conductance parameters underlying 
the Ym fluctuations of the up-states of the slow oscillation, as well as those underlying the continuous 
fluctuations following PPT stimulation. In both cases, the average estimated inhibitory conductance gjo 
was markedly higher than the average estimated excitatory conductance ^^0^ and similarly the estimated 
variance of inhibition o? was higher than the variance of excitation Og. 

The second study shows similar results across the natural wake-sleep cycle of the cat (Fig. [2]): for a major- 
ity of cells, especially during slow-wave-sleep up-states, inhibition dominated in terms of both mean and 
variance. At the population level, the ratio of inhibition to excitation was higher during slow-wave-sleep 
up-states compared to the wake state. In 3 neurons that were recorded across several states, both aver- 
age conductances together with their variances decreased in the wake state compared to slow-wave-sleep 
up-states. In addition, especially during the wake state, some cells displayed comparable excitation and 
inhibition or even a dominant excitation (2 out of 11 cells in the wake state). The study also reports an 
important diversity in the absolute values of the estimated conductance parameters. 

An important concern in this type of studies is the estimation of the leak parameters of the recorded neurons. 
The down-states of the slow oscillation, when the local network is presumably silent, are too short for 
properly estimating these parameters, since these brief periods immediately following prolonged up-states 
are likely to include after-hyperpolarizing currents (Sanchez- Vives et al., 2000) that would bias the estimate 
of the leak. In both cases, values obtained from previous work (Pare et al, 1998; Destexhe and Pare, 1999) 
were used for the up-states: these studies evaluated that the ratio of input resistance after TTX block 
of synaptic activity to the total input resistance during up-states was about 4 to 6-fold. The total input 
resistance was estimated in all cells studied from the linear portion of the I-V curve obtained during up- 
states. During wake or wake-like states, the mean total input resistance, as compared to the mean total input 
resistance during up-states, was taken into account to predict the ratio of input resistance with and without 
synaptic activity. The underlying hypothesis for these assumptions is that the ratio of leak conductance 
to synaptic conductance is similar in all cells during comparable network states, implying that the leak 
conductance and the total synaptic conductance covary in a strong way across cells. In the natural wake- 
sleep cycle study, the dependence of the conductance estimates on this important ratio was systematically 
evaluated: this analysis showed that the estimated domination of inhibition was qualitatively robust, and 
that this prediction failed only when the leak conductance and the total synaptic conductance were assumed 
to be approximately equal. 

On the other hand, the variance estimates do not depend on the leak conductance, provided that the total 
conductance is known, as is the case in the cited studies. They are, however, dependent on the membrane 
capacitance C (see Eqs.|9]and 10). This value was assumed to be constant across the cells analyzed in vivo 



(supposing a specific membrane capacitance of 1 /jp/crn^ and a membrane area of around 30000 /jm^). In 
future studies, it could possibly be estimated using short current pulse injection during the spontaneous 
activity, where the capacitance can be extracted from the time constant of exponential fits to the V,„ decay. 



Synaptic conductance variances are useful parameters that can be related to spike initiation (see Sections 3.3 



and 3.4). The VmD method provided the first published estimates for these parameters in vivo (for a 



different approach, see Monier et al., 2008, this issue). 
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3.2.1 Outline of the method 

The point-conductance model given by Eqs. [l]-[3] was studied further, and recently we showed that the 
power spectral density (PSD) of the fluctuations described by this model can be well approximated by 
the following expression (Destexhe and Rudolph, 2004): 



4 



1 +0)2x2 



+ 



1 + 0)2 x2 1 + 0)2 x2 



(11) 



where (0 = 2nf, f is the frequency, Gj = Gl + geO+ gio is the total membrane conductance, v = C/Gj is 
the effective time constant, and V = {GlEl + geoEe + gioEi) /Gt is the average membrane potential. The 
"effective leak" approximation used to derive this equation consisted in incorporating the average synaptic 
conductances into the total leak conductance, and then considering that fluctuations around the obtained 
mean voltage are subjected to a constant driving force (Destexhe and Rudolph, 2004). 

As mentioned above, the synaptic time constant parameters, x^, and x,-, need to be estimated in order to ex- 
tract precise values of the conductance variances with the VmD method. As those two parameters appear in 
the theoretical expression of the PSD, we explored the possibility of evaluating them from the analysis 
of the PSD of experimentally recorded Vm fluctuations. To this end, the following simplified expression 
can be fitted: 



1 + 0)2 x2 



Ag Xg Ai X,' 



+ 



(12) 



1 + 0)2x2 l+a)2x2_ 

where Ag and A, are amplitude parameters. This five parameter template is used to provide estimates of the 
parameters Xg and x, (supposing that x„, has been measured). A further simplification consists in assuming 
that Ae = Aj, which was used for fitting in vivo data (see Methods). 

These analytic expressions were tested by comparing the prediction to numerical simulations of a single- 
compartment model subject to fluctuating synaptic conductances (Eqs. [l]-[3]). The matching between the 
analytic expression and the PSD obtained numerically was nearly perfect, as shown in Fig. |3]^ and as 
reported in detail previously (Destexhe and Rudolph, 2004). 



3.2.2 Testing synaptic time constants estimates with dynamic-clamp 

We applied the procedure described above to the PSD of Y,„ fluctuations obtained by controlled, dynamic- 
clamp fluctuating conductance injection in cortical neurons in vitro (using a new, high resolution electrode 
compensation technique, Brette et al., 2005, 2007; Rudolph et al., 2005). In this case, the scaling of the 



PSD conforms to the prediction (Fig.pB): the theoretical template (Eq. 12) can provide a very good fit of 



the experimentally obtained PSD, up to around 400 Hz, where recording noise becomes important. The 



template used was according to Eq. [TT]or Eq. [T2[ both of which provided equally good fits (not shown). 
This shows that the analytic expression for the PSD is consistent not only with models, but also with 
conductance injection in real neurons in vitro. 
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We have also attempted to apply the same procedure to V„, fluctuations resulting from real synaptic activ- 
ity, during up-states recorded in vitro (Fig.[3p) and during sustained network activity in vivo (Fig.[3p). In 
this case, however, it is apparent that the experimental PSDs cannot be fitted with the theoretical template 
as nicely as for dynamic-clamp data (Fig. [3^). The PSD presents a frequency scaling region at high fre- 
quencies, and scales as 1//" with a different exponent a as predicted by the theory (see Figs.|3p-D). The 
analytic expression (Eq. 1 1 1 predicts that the PSD should scale as 1/f^ at high frequencies, but the exper- 
iments show that the exponent a is obviously lower than that value (see Discussion for possible reasons 
for such a difference). This difference of course compromises the accuracy of the method to estimate Xg 
and X, in situations of real synaptic bombardment. Nevertheless, including the values of x^ = 3 ms and 
X, = 10 ms provided acceptable fits to the low-frequency (<100 Hz) part of the spectrum (Fig. |3]C-D, red 
curves). However, in this case, small variations (around 20-30%) around these values of Xg and X/ yielded 
equally good fits (not shown; see also Rudolph et al., 2005). Thus, the method cannot be used to precisely 
estimate those parameters, but can nevertheless be used to broadly estimate them with an error of the order 
of 30%. 



3.3 Estimating spike-triggering conductance configurations 
3.3.1 A preliminary investigation 

The conductance measurements outlined above show that there is a diversity of combinations of ge and gi 
that underlies the genesis of subthreshold activity in different preparations. We used a computational model 
based on the point-conductance model, but including a spiking mechanism (see Eq.|4]), to reproduce those 
measurements and try to infer what different properties such states may have in terms of spike selectivity. 
Indeed, an infinite number of combinations of ge and gj can give similar V„, activity. Figure |4j\ illustrates 
two extreme examples out of this continuum: first a state where both excitatory and inhibitory conductances 
are of comparable magnitude (Fig. |4j\, left; "Equal conductances"). In this state, both conductances are 
lower than the resting conductance of the cell and the is fluctuating around -60 mV. Second, similar 
fluctuations can be obtained when both conductances are of larger magnitude, but in this case, inhibition has 
to be augmented several-fold to maintain the around -60 mV (Fig.|4j\, right; "Inhibition-dominated"). 
Such conductance values are more typical of what is usually measured in vivo (Rudolph et al., 2005, 2007; 
Monier et al., 2008, this issue). Both conductances are larger than the resting conductance, a situation 
which can be described as a "high-conductance state". 

To determine how these two states differ in their spike selectivity, we evaluated the spike-triggering con- 
ductances by averaging the conductance traces collected in 50 ms windows preceding spikes. This average 
pattern of conductance variations leading to spikes is shown in Fig.|4^. For equal-conductance states, there 
is an increase of total conductance preceding spikes (purple curve in Fig.|4^, left), as can be expected from 
the fact that excitation increases (ge curve in Fig.|4^, left). In contrast, for inhibition-dominated states, the 
total conductance decreases prior to the spike (purple curve in Fig.|4^, right), and this decrease necessarily 
comes from a similar decrease of inhibitory conductance, which is, in this case, stronger than the increase 
of excitatory conductance (gi curve in Fig.|4^, right). Thus, in such states the spike seems primarily caused 
by a drop of inhibition. 

This pattern was seen not only in the average, but also at the level of single spikes. Using a vector represen- 
tation to display the conductance variation preceding spikes (each vector links the conductance state in a 
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window of 30-40 ms before the spike with that in the 10 ms preceding the spike) shows that the majority of 
spikes follow the average pattern (Fig.|4p). The same features were also present when the integrate-and-fire 
model was used (not shown), and thus do not seem to depend on the spike generating mechanisms. 

These patterns of conductance variations preceding spikes were also investigated in real neurons by using 
dynamic-clamp experiments to inject fluctuating conductances in vitro. In this case, performing the same 
analysis as above revealed similar features: spike-triggered averages (STAs) of the injected conductances 
displayed either increase or decrease in total conductance, depending on the conductance parameters used 
(Fig.[5}\), and the vector representations were also similar (Fig. [5^). It suggested that these features are 
independent of the spike generating mechanism but rather are caused by subthreshold dynamics. 



3.3.2 A geometrical interpretation based on the point-conductance model 

The configuration of synaptic conductances just before spikes can be explained qualitatively by considering 
that the total current must be positive at spike time, i.e., geiEg — Vt) +gi{Ei — Vt) + Gl{El — Vt) > 0, where 
Vt is the spike threshold (using an integrate-and-fire approximation). 

This inequality defines a half-plane in which {ge-,gi) must lie at spike time. Fig. [6|^. shows graphically 
how this inequality affects the synaptic conductances. The variable (ge.gi) is normally distributed, so that 
isoprobability curves are ellipses in the plane (plotted in red). In that plane, the line {ge + gi = geO + 
gio} going through the center of the ellipses defines the points for which the total conductance equals 
the mean conductance, and the line {ge{Ee — Vt) +gi{Ei — Vt) + Gl{Ei — Vt) = 0} defines the border of 
the half-plane in which conductances lie at spike time. In the equal conductances regime (Fig. |6]\, left), 
synaptic conductances are small and have similar variances, so that isoprobability curves are circular; the 
intersection of the half-plane with those circles is mostly above the mean total conductance line, so that the 
total conductance is higher than average at spike time. 

In the inhibition-dominated regime (Fig.[6j^, right), synaptic conductances are large and the variance of gi is 
larger than the variance of ge, so that isoprobability curves are vertically elongated ellipses; the intersection 
of the half-plane with those ellipses is essentially below the mean total conductance line, so that the total 
conductance is lower than average at spike time. 

More precisely, when isoprobability curves are circular (equal variances), then the expected total con- 
ductance is unchanged at spike time when the lines {ge{Ee — Vt) + gi{Ej — Vt) + Gi{Ei — Vt) =0} and 
{ge + gi = geO + gio} ^rc orthogonal, i.e., when Ee — Vt + Et — Vt = 0. Spikes are associated with increases 
in conductance when the first line has a higher slope, i.e., when Ee — Vt> Vt — Et (which is typically the 
case). 

When isoprobability curves are not circular, we can look at the graph in the space {^,§:) where isoproba- 
bility curves are circular. Then the orthogonality condition between the lines 

( ^OeiEe - Vt) + ^o,{Ei - Vt) + Gl(El - Vt) = o] 

and 

/^^^ , "^'^ , \ 

< Oe H Oi = geO+giO (■ 

[ <3e <3i } 

reads 

ol{Ee-Vt)+oUE-Vt)=0. 
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It follows that spikes are associated with increases in total conductance when the following condition is 
met: 

\Ie,-V,- 

One can also recover this result by calculating the expectation of the conductance change conditionally 
to the fact that the current at spike threshold is positive (implicitly, we are neglecting the correlation time 
constants of the synaptic conductances). Using typical values (Vt = —55 mV, Eg = mV, Ei = —75 mV), we 
conclude that spikes are associated with increases in total conductance when Og > 0.6o,-. This inequality is 
indeed satisfied in the equal conductances regime and not in the inhibition-dominated regime investigated 
above. 

3.3.3 Testing the geometrical prediction with dynamic-clamp 

The geometrical reasoning predicts that the sign of the total synaptic conductance change triggering spikes 
depends only on the ratio of synaptic variances, and not on the average conductances. We have systemati- 
cally tested this prediction using dynamic-clamp injection of fluctuating conductances in vitro. In 8 regular 
spiking cortical neurons, we scanned different parameter regimes in a total of 36 fluctuating conductance 
injections. Fig.|6^ shows two examples from the same cell: both correspond to an average "high conduc- 
tance" regime, dominated by inhibition, but in one case it is the variance of excitation, in the other case 
the variance of inhibition, that is higher. We can see that the total conductance before the spike increases 
in the first case, but decreases in the second. Fig. |6p (left) shows the average total conductance change 
preceding spikes as a function of Og/o,, for all the 36 injections: the vertical dashed line represents the 
predicted value of Oe/o,=0.6, which indeed separates all the "conductance drop" configurations from the 
"conductance increase" configurations. Even though the prediction is based on a simple integrate-and-fire 
extension of the point-conductance model, we can see that the ratio of synaptic variances can predict the 
sign of the total conductance change triggering spikes in biological cortical neurons subjected to fluctuating 
excitatory and inhibitory conductances. 

In addition (Fig. |6p, right), the dynamic-clamp data shows that the average amplitude of change (Age = 
geo kg, see Methods) of each synaptic conductance preceding a spike is related, in a linear way, to the 
standard deviation of this conductance. For a fixed value of standard deviation, there was no significant 
influence of the average conductance (not shown). This observation is consistent with the idea that in all 
the cases studied here, the firing of the cell was driven by fluctuations in the V,„, rather than by a high mean 
Wfn value (not shown). 

Taken together, these theoretical and experimental analyses indicate that the average total conductance drop 
preceding spikes, as seen in the "high conductance" case we initially considered (Fig. |4]), is not a direct 
consequence of the "high conductance" state of the membrane, but is in fact related to the high inhibitory 
variance, which is indeed to be expected especially when the mean inhibitory conductance is also high (as 
confirmed by the studies presented in the first part of this article). 

3.4 Estimating spike-triggered averages of synaptic conductances from the \rn 

In order to extend the analysis of spike-triggering conductance configurations to real recurrent cortical 
activity observed in vivo, we recently developed a procedure to extract the spike-triggered averages (STAs) 
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of conductances from recordings of the Y,n (Fig.|7| Pospischil et al., 2007). The STA of the V,,, is calculated 
first, and the method searches for the "most likely" spike -related conductance time courses (geit), gi{t)) 
that are compatible with the observed voltage STA. The procedure is based on a discretization of the time 
axis in Eqs.[l]-[3| which, rearranged, lead to the following relations: 



V'-Ei 1 XL C At C 



= (15) 



o,- V 2A? V ' ^ / X, 



where Xl = C / Gl is the resting membrane time constant. Note that and have become Gaussian- 
distributed random numbers and . There is a continuum of combinations {ge^^, g^^^} that can advance 
the membrane potential from y^+^ to y^+^, each pair occurring with a probability 

/ := Pter',sr'lsf,sf) = 4^-5e'+«'' = ^.-=''', (16) 



Because of Eq. 13. and g'^ are not independent and is, thus, a unidimensional distribution only. 



Given initial conductances {gg,^^}, one can write down the probability p for certain series of conductances 
{ge.gj} j=o,...,n to occur that reproduce a given voltage trace {y'}/=i....,;i+i: 

n-l 

P=YIp' (18) 

Due to the symmetry of the distribution p, the average paths of the conductances coincide with the most 
likely ones. It is thus sufficient to determine the conductance series with extremal likelihood by solving the 
n-dimensional system of linear equations 



0| , (19) 

e ) k=\,...,n 



dge 



where X = Y!k=o^^' ^'^^ vector {gg}. This is equivalent to solving = 0}/t=i,...,n and involves the 
numerical inversion of an n x n-matrix, which can be done using standard numeric methods (Press et al., 
1986). The series {g^} is subsequently obtained from Eq. 13 Details of this procedure as well as an 



evaluation of its performance can be found in Pospischil et al., 2007. 

The method requires first an estimation of the parameters describing the distribution of each of the con- 
ductances, which can be obtained by the VmD method. The leak parameters of the cell, or alternatively 
the effective parameters during V„, fluctuations (effective conductance and effective time constant), also 
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have to be estimated prior to this analysis. The method was successfully tested using computational mod- 
els: conductance STAs extracted from the Y,n of an integrate-and-fire point-conductance model are nearly 
identical to the numerically obtained conductance STAs (Fig.|7]A, left; Pospischil et al., 2007). 

During response to sensory stimuli, there can be a substantial degree of correlation between excitatory and 
inhibitory synaptic input (Monier et al., 2003; Wehr and Zador, 2003; Wilent and Contreras, 2005). Since 
this situation has not been addressed in Pospischil et al., 2007, we would like to sketch a possible extension 
of the method. To this end, we reformulate the discretized versions of Eqs.|2||3]in the following way: 

k+l _ k 
At 

At 



?e-SeO 



?'i-giO 



2At 
x,(l+c) 

2At 



&i + V~c^2). 



(20) 
(21) 



Here, instead of having one "private" white noise source feeding each conductance channel, now the same 
two noise sources and ^2 contribute to both inhibition and excitation. The amount of correlation is tuned 
by the parameter c. Also, since there is evidence that the peak of the ^e-^,-crosscorrelation is not always 
centered at during stimulus-evoked responses ("delayed inhibition"; see Wehr and Zador, 2003; Wilent 
and Contreras, 2005), we allow a non-zero delay d: for a positive parameter d, the inhibitory channel 



receives the input that the excitatory channel received d time steps before. Eqs. 20 and 21 can be solved for 
and thus replacing Eqs. 14 andflsl It is then possible to proceed as in the uncorrelated case, where 



now, due to the delay, the matrix describing Eq.[T9]has additional subdiagonal entries. 

However, the application of this extended method requires the estimation of the usual leak parameters, of 
conductance distribution parameters - for which the VmD method cannot be directly used in its current 
form since it is based on uncorrelated noise sources - as well as knowledge of the parameters c and d. At 
present, we can only speculate on how c and d could be evaluated in experiments: extracellularly recorded 
spike trains could perhaps be used to this end, provided that simultaneously recorded single units could be 
classified as excitatory or inhibitory. Alternatively, different plausible c and d values could be scanned to 
examine how they could potentially influence the conductance STAs extracted from a given STA. 



3.4.1 Testing STA estimation with dynamic-clamp 

The dynamic-clamp data presented above was used to evaluate the accuracy of the conductance STA es- 
timation method (for the uncorrelated case only): indeed, the conductance STAs estimated from the Vm 
STAs could be compared to conductance STAs obtained directly by averaging conductance traces, since 
in dynamic-clamp the injected conductances are perfectly controlled by the experimentalist. For the 36 
injections analyzed, a good match was observed between the two (see one example in Fig.|7}A, right panel). 
Since we intended to evaluate the STA method specifically, we assumed that conductance distribution pa- 
rameters were known. In addition, we estimated the cell's leak parameters (see Methods). To quantify 
the comparison on a population basis, we did the following analyses: exponential functions (see Methods) 
were fitted to each conductance STA (estimated and directly measured) starting at 1 ms before the spike 
and decaying to baseline backwards in time. We then compared the asymptotic values (i.e., the average 
baseline conductances) and the time constants of these fits, as well as the amplitude of conductance change 
from the start of the fit to the asymptote (i.e., the amplitude of conductance change preceding the spike). In 
all cases but one, the exponential functions provided excellent fits to the conductance STAs. It was neces- 



Piwkowska et ah, J. Neurosci. Meth. (in press, 2008) 



16 



sary to exclude the 1 ms time window preceding the spike to avoid severe contamination of the analyses by 
intrinsic conductances. Excluding a broader time window did not improve the analyses for these neurons. 

The average baseline conductances always matched very well (error of -0.1 ± 0.25 nS, or -0.8 ± 2.6 %, 
for excitation; 0.8 ±1.5 nS, or 0.6 ± 4.5%, for inhibition; not shown). More importantly, the estimates of 
the average conductance patterns leading to spikes were also in good correspondence with the measured 
patterns, in terms of both the amplitude of conductance change (Fig.|7p, top; error of -1.2 ± 3 nS, or -26 
± 28.8 %, for excitatory amplitude change, and -2.0 ± 2.5 nS, or -10.7 ± 47%, for inhibitory amplitude 
change) and the time constant (Fig.|7]C, bottom; error of 0.39 ± 0.48 ms, or 11.2 ± 21.1% for excitatory 
time constant, and 0.36 ± 1.71, or 2.6 ± 18.8%, for inhibitory time constant). For excitation, the error on 
the estimate of the amplitude is correlated with the error on the estimate of the time constant (not shown): 
this suggests that, in most cases, a slightly too fast rise of the excitatory conductance results in a slightly too 
high amplitude of conductance change. The lack of correlation between the two error measures in the case 
of the inhibitory conductance points to a more complex origin for the observed errors. Moreover, the errors 
on the amplitude of the two conductance changes are positively correlated (Fig.|7p), top). This dependency 
actually ensures that the error on the estimated total conductance change (excitation-inhibition) remains 
small (-0.8 ± 2.4 nS; Fig.|7p, bottom). 

Finally, we have investigated the dependency of the estimate errors on a diversity of variables, and found 
a correlation of amplitude errors with the average V,,, during the fluctuating conductance injection (not 
shown): this dependency points to a possible contamination by intrinsic conductances activated differen- 
tially at different average levels. However, we found no dependency of the error on average firing rate 
(for the rates up to around 30 Hz studied here), suggesting no important contamination by spike-dependent 
conductances like the ones underlying the after-hyperpolarization (when care is taken to compute STAs 
using spikes preceded by at least 100 ms of silence, see Methods). 

We have also verified that conductance STA estimates can be relied on to investigate what factors determine 
the average conductance variations preceding spikes. Fig. [8] shows that the analyses performed previously 
on dynamic-clamp data (i.e., on STAs obtained directly by averaging the conductance traces. Fig. [6^-C), 
can also be successfully performed using the conductance STAs estimated from the corresponding 
STAs: Fig. [8]^ shows the estimated STAs in the two different "high conductance" states, dominated by 
either excitatory or inhibitory variance (compare to Fig. [6p). Fig. [8p shows, for the population data, the 
(significant) correlations between total conductance change and Og/a,, as well as between the change of 
each of the conductances and the corresponding standard deviation. Note the similarity between Fig. [6^-C 
and Fig. [8|\-B, even though the correlations at the population level are more noisy when the estimated 
STAs are used. 

3.4.2 STA Analysis of intracellular recordings of cortical neurons in vivo 

The conductance STA estimation method was used to determine conductance variations preceding spikes 
during fluctuations in vivo (Rudolph et al., 2007). Starting from recordings of spontaneous spiking 
activity in awake or naturally sleeping cats, we computed the spike-triggered average of the Y,„ (Fig. |9]). 
Using values of geo, gio, Og, o, estimated using the VmD method (see above), we computed the most likely 
conductance traces yielding the observed Y,n averages. Most of these analyses (7 out of 10 cells for awake, 
6 out of 6 for slow-wave-sleep, 2 out of 2 for REM) revealed conductance dynamics consistent with states 
dominated by inhibitory variance: there was a drop of the total conductance preceding spikes, due to a 
strong decrease of the inhibitory conductance (Fig.|9} right). However, a few cases, in the wake state (3 out 
of 10 cells), displayed the opposite configuration with the total synaptic conductance increasing before the 
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spike (Fig. |9| left). 

We also checked how the geometrical prediction relating the sign of total conductance change preceding 
spikes and the ratio Oe/a, performed for this data (Fig. |9^). We have seen that the critical value of Og/a/ for 
which the total conductance change shifts from positive to negative depends on the spike threshold. This 
parameter was quite variable in the recorded cells, and so a critical Oe/O/ value was calculated for each cell. 
Fig.|9p shows the lowest and highest critical values obtained (dashed lines), and also displays in white the 
cells which do not conform to the prediction based on their critical value. This is the case for only 4 out of 
18 cells, for three of which the total conductance change is close to zero. 

The extraction of conductance STAs depends on the accuracy of the synaptic conductance parameters 
estimated with the VmD method, which means that assumptions made about the leak conductance and the 



cell capacitance will influence the results (see Section 3.1.3). We have shown that the ratio Oe/a, should 
determine whether, on average, the total conductance increases or decreases prior to the spike. This ratio is 
independent of the leak conductance, but it depends on the capacitance C. However, it appears in both the 
numerator and the denominator of the ratio. Fig. |9p shows the dependency of Og/o, on this parameter for 
different values of total input resistance and two sets of realistic values for the V,,, distribution parameters. 
This analysis indicates that a reasonable error on C produces a limited error on the ratio Oe/o,, and suggests 
that conclusions drawn from the in vivo data about the respective contributions of excitation and inhibition 
in triggering spikes are valid. 



4 Discussion 



4.1 Synaptic conductance analysis methods based on the fluctuating point-con- 
ductance model 

We presented how the simple point-conductance model of cortical synaptic activity can provide a basis for 
the analysis of experimental data, essentially through the matching of expressions derived from the model 
to intracellular recordings of cortical neurons. This approach has been used for extracting different 
parameters from the recurrent cortical activity in vivo: the averages and variances of excitatory and in- 
hibitory conductances, their decay time constants and the optimal conductance waveform underlying spike 
selectivity. These analyses were possible because the point conductance model represents in a compact and 
mathematically tractable way the activity resulting from several thousand synapses. 

The VmD analysis provides a characterization of synaptic activity in simple terms (average conductance, 
level of fluctuations). Such parameters can readily be incorporated in computational models to yield the 
Vm and conductance state corresponding to in vivo activity with just a few variables. This approach has 
been used for example in network simulations to obtain realistic conductance states in neurons even with 
small networks (Haeusler and Maass, 2007). It is also directly usable in dynamic-clamp experiments to 
investigate the impact of synaptic background activity on signal processing by single cortical or thalamic 
neurons (Fellous et al., 2003; Shu et al., 2003; Wolfart et al., 2005; Desai and Walcott, 2006). 

Beyond the matching of experimental distributions to a theoretical expression (VmD method), we 
have also attempted to match the PSDs of V^ fluctuations. This approach provides some validation for 
assumptions made about synaptic time constants on the basis of published studies (Destexhe and Pare, 
1999; Destexhe et al., 2001). However, the fact that the point-conductance model does not account for 
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the scaling properties of experimental PSDs (Fig.[3p-D) limits the accuracy of the method and yields only 
broad estimates of the synaptic time constants (approximately 30% error). A parallel study (Bedard and 
Destexhe, 2007) has shown that the frequency scaling observed experimentally cannot be accounted for by 
standard cable theory, and modifications of cable equations are required to match those values. We hope to 
obtain a more accurate fitting template, which would allow more precise PSD analyses in the future. 

Finally, in order to be able to study average spike-triggered patterns of conductances in vivo, we have re- 
cently developed a probabilistic method for extracting the STAs of conductances from STAs of the 
(Pospischil et al., 2007). This method relies on the estimation of a number of parameters inherent to the 
neuron (like leak conductance and capacitance), as well as on synaptic conductance parameters estimated, 
for example, with the VmD method. We provided new data from dynamic-clamp experiments in vitro, 
demonstrating that given the synaptic conductance parameters, a good match between the extracted con- 
ductance STAs and the actual, known conductance STAs can be obtained. 

4.2 Comparison to other synaptic conductance analysis methods 

Other studies proposed synaptic conductance estimates in vivo derived from I-V curves (in current-clamp) 
or V-I curves (in voltage-clamp) obtained at different points in time following a stimulus (Borg-Graham et 
al.,1998; Anderson et al., 2000; Monier et al., 2003; Wehr and Zador, 2003; Wilent and Contreras, 2005) 
or the onset of an up-state (Haider et al., 2006). An exhaustive comparison of these studies is beyond the 
scope of the present article (see Monier et al., 2008, this issue), but a few points can be stressed. The fact 
that the VmD method applies to current-clamp data obtained at a few levels of constant injected current 
circumvents technical problems with voltage-clamp due to high series resistance of patch electrodes in 
vivo (but see Borg-Graham et al., 1998; Monier et al., 2003; Wehr and Zador, 2005) or the need for 
discontinuous voltage-clamp with sharp electrodes (Haider et al., 2006). The fact that it relies on a strong 
assumption about the stochasticity of synaptic inputs and the independence of excitation and inhibition 
allows the analysis of spontaneous cortical activity with no "zero" time point. However, this assumption 
also makes it unsuited, in its present form, for the analysis of stimulus-evoked activity with important 
temporal structure. 

The VmD method suffers from one common limitation with other approaches for synaptic conductance es- 
timation: the need to separate synaptic currents from leak currents. The fact that it is used for the analysis 
of on-going, spontaneous activity in cortical networks poses, however, an additional complication: indeed, 
when synaptic inputs evoked by sensory stimulation are analyzed, it is with reference to the pre-stimulus 
activity, which includes the leak current and any other baseline currents, including on-going synaptic ac- 
tivity (Borg-Graham et al.,1998; Anderson et al., 2000; Monier et al., 2003; Wehr and Zador, 2003; Wilent 
and Contreras, 2005). When we attempt to analyze the spontaneous activity itself, there is no straightfor- 
ward reference that can be used. As mentioned above, the short down-states do not seem a good candidate 
since they include after-hyperpolarizing currents consecutive to the up-states (Sanchez- Vives et al., 2000), 
although a down- state-referenced analysis (as in Haider et al., 2006) could perhaps be compared to the 
analysis performed in the reviewed studies. However, in the wake state, the continuous on-going activity 
does not even present down- states. The only precise approach to evaluate the leak conductance for each 
studied cell is to block all synaptic activity with TTX, but that also means recording only one or two cells 
per animal for this protocol, which is an extremely constraining experimental situation. The compromise 
chosen in the reviewed in vivo studies consisted in using published average values obtained in previous 
TTX experiments in a similar preparation (Pare et al., 1998) and checking the robustness of the estimates 
to the assumed leak conductance parameter. The estimates should be re-evaluated in the light of any future 
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experimental data providing information about the leak parameters of cortical cells in vivo in the absence 
of synaptic activity. Optimally, such future studies would include the simple, classical protocol required 
for a subsequent VmD analysis: several-second- long current-clamp recordings of the at different levels 
of steady injected current. 

As to our probabilistic method for extracting the STAs of conductances from STAs of the Vm (Pospischil et 
al., 2007), we are not aware of any other method currently allowing this analysis for spontaneous activity: in 
voltage-clamp, no spikes are recorded and so, obviously, conductances leading to spikes cannot be extracted 
directly. Voltage-clamp can only be used for extracting plausible conductance STAs when the precise times 
of spikes are known and reproducible (for example, at a given delay after a sensory stimulation), so that 
the estimated, stimulus-locked conductance dynamics can be reasonably expected to lead to spikes in the 
current-clamp configuration (Monier et al., 2003; Wehr and Zador, 2003). 

4.3 Dependency of the STA extraction method on the VmD method 

The result of the STA analysis is, however, dependent on the estimates of synaptic averages and variances, 
so that we may ask to what extent this result depends on the accuracy of the VmD method. In general, 
the average baseline conductances will reflect the estimates of average synaptic conductances (since both 
analyses are constrained by the same total input resistance measure and the same leak conductance as- 
sumption), unless the total input resistance changes markedly between the current levels used for the VmD 
analysis and the zero-current level at which STAs are extracted, due to activation of intrinsic conductances: 
in this case, since the conductances (Gl, geit), gii^)) have to be compatible with the voltage V„i{t), the 
conductance STA baselines can be considerably shifted away from the mean conductance values tens of ms 
before the spike. The same effect is seen in in vitro dynamic-clamp experiments when a wrong leak con- 
ductance value is used (not shown). This distortion could be used as an indication for important activation 
of intrinsic conductances and suggest that the result should be discarded. 

The variances of the synaptic conductances are correlated to the amplitudes of average conductance change 
preceding a spike (as shown in Fig. [6]). The estimation of synaptic conductance variances with the VmD 
method is independent of the assumption made about the leak conductance, which excludes this source 
of potential error. It is dependent on the membrane capacitance C, which should be evaluated on a cell- 
by-cell basis whenever possible in future studies. However, we have also shown that the ratio of synaptic 
conductance variances is only weakly dependent on the precise value of C: this ratio determines the sign of 
total conductance change preceding a spike, so that the estimates of this sign from in vivo data seem robust 
(see below). 

4.4 Dynamic-clamp as a tool to evaluate conductance analysis methods 

We illustrated how and to what extent the validity of different approaches for conductance analysis can be 
tested using dynamic -clamp. This electrophysiological technique is an attractive tool to evaluate methods of 
conductance analysis, since it allows to mimic the activation of known conductances in a biological neuron: 
the results of an analysis method based on V,„ recordings can be directly compared to measures of the actual 
conductances controlled by the experimentalist. In all the dynamic-clamp applications presented here, we 
have used the same description for the synaptic conductances - the Ornstein-Uhlenbeck stochastic model - 
as in the theoretical analyses. This means that we could compare how the analysis methods perform if the 
stochastic conductances are inserted at the soma of a real cortical neuron, with a complex structure and a 
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variety of intrinsic channels, instead of a passive single compartment. As we have seen, potential dendritic 
effects solicited only during distributed synaptic stimulation, like, possibly, the unexpected scaling of the 
V,„ PSDs during real synaptic activity (Bedard and Destexhe, 2007), cannot be addressed with this somatic 
injection technique. They could perhaps be investigated in the future using dendritic patch-clamp. 

The comparison was performed most extensively for the conductance STA estimation method: it indicates 
that if the window of analysis is chosen properly (by excluding a window of about 1 ms before the spike, 
and also excluding inter-spike-intervals shorter than around 100 ms), the estimations perform well, and 
that the estimation errors, correlated with the average V,„, are presumably linked to V,„-dependent intrinsic 
channels. However, we did not systematically compare how the methods perform if the synaptic conduc- 
tances deviate from the Omstein-Uhlenbeck model: such an approach could constitute another application 
of the dynamic-clamp tool to the evaluation of conductance analysis methods. We have also not attempted 
yet to evaluate in dynamic-clamp the extended STA analysis method sketched above, which incorporates a 



known correlation between excitation and inhibition (Eqs. 20-21 1. 



4.5 Patterns of excitation and inhibition triggering spikes depend on the variances 
of the synaptic inputs 

What are the patterns of excitation and inhibition triggering spikes under different conditions of network 
activity ? A preliminary study, using both models and dynamic-clamp conductance injection in vitro, 
pointed to the fact that these patterns depend on the statistics of synaptic conductances, and that spikes could 
be preceded, on average, either by increases in total synaptic conductance, indicating a predominant role for 
excitation, or by decreases in total synaptic conductance, indicating a predominant role for inhibition. Other 
authors (Hasenstaub et al., 2005) have recently suggested that drop of inhibition can play an important role 
in determining spike timing in cortical neurons, based on dynamic-clamp injection of specific synaptic 
conductances with parameters matched to their in vivo recordings. Here we explored this issue further 
and showed, first by a theoretical reasoning, and second by scanning different parameter regimes using 
dynamic-clamp conductance injection, that the sign of the total conductance change before a spike does 
not directly depend on the average synaptic conductances, but is solely determined by the ratio of synaptic 
conductance variances. The variance of synaptic conductances in the Ornstein-Uhlenbeck model is related 
to the degree of correlation between Poisson input trains of each type in a detailed biophysical model 
(Destexhe et al., 2001): the level of synchrony among inhibitory neurons could thus be determining for 
spike timing whenever it significantly exceeds the level of synchrony among excitatory neurons. The 
rule relating the ratio of synaptic conductance variances with the sign of the average conductance change 
preceding a spike appears as a generalization of more specific results obtained by Tateno and Robinson 
in cortical neurons (Tateno and Robinson, 2006), as well as by the group of Jaeger in the cerebellum 
(Gauck and Jaeger, 2000, 2003; Suter and Jaeger, 2004): these authors also used dynamic-clamp injection 
of synaptic conductances and the conductance STAs they obtained seem consistent with our rule. 

This result stresses the importance of evaluating the variances of synaptic conductances, in addition to 
their averages, when analyzing Ym fluctuations recorded in vivo. To our knowledge, the two reviewed 
studies using the VmD method (Rudolph et al., 2005, 2007) are the only ones (together with Monier et al., 
2008, this issue) explicitely providing estimates of synaptic conductance variances. The application of the 
probabilistic method also allowed for the first time the extraction of synaptic conductance STAs from in 
vivo recordings of spontaneous fluctuations in awake and naturally sleeping cats (Rudolph et al., 2007), 
and led to the observation of both types of firing regimes described above - average total conductance 
increase and average total conductance drop - with a majority of cases displaying the inhibition-dominated. 
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Figure 1: Conductance extraction from V,„ distributions: numerical and dynamic-clamp test of the method. 
A. Simulation of the point-conductance model (top trace) and comparison between numerically computed 
Vm distributions (bottom; blue) and the analytic expression (red curves; conductance values shown in 
the bar graph). B. Dynamic-clamp injection of the point-conductance model in a real neuron. (Right) 
Conductance parameters are re-estimated (back, colored; error bars are standard deviations obtained when 
the same injected conductance parameters are re-estimated in different cells) from the V„j distributions and 
compared to the known parameters of the injected conductances (front, grey). (Left) The experimental 
Y,n distributions are compared to the analytic distributions calculated using the re-estimated conductance 
parameters. C. Comparison of a spontaneous up-state (Natural up-state) with an artificial up-state recreated 
using conductance injection (Dynamic-clamp). Modified from Rudolph et al., 2004. 
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Figure 2: VmD estimation of conductances from intracellular recordings in awake and naturally sleep- 
ing cats. A. Intracellular recordings in awake and naturally sleeping (SWS) cats. Recordings were made 
in association cortex (area 5-7). B. Examples of V,„ distributions computed during wakefulness (Awake) 
and slow-wave sleep up-states (SWS). The continuous lines show Gaussian fits of the experimental dis- 
tributions. Insets: current-voltage relations obtained for these particular neurons. C. Conductance values 
estimated using the VmD method. Results for the means (geo, gjo) and standard deviations (Oe, o,) of 
excitatory and inhibitory conductances, respectively, as well as their ratios are shown (error bars: standard 
deviations obtained by repeating the analysis using different pairs of injected current levels). D. Grouped 
data showing the means and standard deviations of the conductances for different cells across different 
behavioral states (REM = Rapid Eye Movement sleep). Figure modified from Rudolph et al., 2007. 
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Figure 3: Fit of the synaptic time constants to the power spectrum of the membrane potential. A. Compar- 
ison between the analytic prediction (Eq. 1 1 , red) and the PSD of the for a single-compartment model 
(Eq. [TJ black) subject to excitatory and inhibitory fluctuating conductances (Eqs. [2]j3| Xg = 3 ms and X/ = 
10 ms). B. PSD of the activity in a guinea-pig visual cortex neuron (black), where the same model 
of fluctuating conductances as in A was injected using dynamic-clamp. The red curve shows the analytic 
prediction using the same parameters as the injected conductances (Xg = 2.7 ms and X/ = 10.5 ms). C. PSD 
of Y,n activity obtained in a ferret visual cortex neuron (black) during spontaneously occurring up-states. 
The PSD was computed by averaging PSDs calculated for each up-state. The red curve shows the best fit 
of the analytic expression with Xg = 3 ms and x, = 10 ms. D. PSD of V„, activity recorded in cat association 
cortex during activated states in vivo. The red curve shows the best fit obtained with Xg = 3 ms and x, = 
10 ms. Panel A modified from Destexhe and Rudolph, 2004; Panel D modified from Rudolph et al., 2005. 
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Figure 4: Comparison between equal conductances and inhibition-dominated states in a computational 
model. A. Equal conductance (left; geo = gio = 10 nS, Og = o, = 2.5 nS) and inhibition-dominated states 
(right; geo = 25 nS, gio = 100 nS, Og = 7 nS and a, = 28 nS) in the point-conductance model. Excitatory 
and inhibitory conductances, and the membrane potential, are shown from top to bottom. Action potentials 
(truncated here) were described by Hodgkin-Huxley type models (Destexhe et al., 2001; Eq.|4]). B. Average 
conductance patterns triggering spikes. Spike-triggered averages (STAs) of excitatory, inhibitory and total 
conductance were computed in a window of 50 ms before the spike. C. Vector representation showing 
the variation of synaptic conductances preceding each spike. The excitatory and inhibitory conductances 
were averaged in two windows of 30-40 ms and 0-10 ms (circle) before the spike, and a vector was drawn 
between the obtained values. 



Piwkowska et ah, J. Neurosci. Meth. (in press, 2008) 



30 



B 



Equal conductances 




Inhibition-dominated 

Total conductance 




40 30 20 10 
Time preceding spike (ms) 



40 30 20 10 

Time preceding spike (ms) 




en 

c 

_o 
15 
c 



40 



30 



20 



10 




2.5 5 7.5 

Excitation (nS) 



5 10 
Excitation (nS) 



Figure 5: Average conductance patterns triggering spikes in dynamic -clamp experiments. A. Spike- 
triggered averages of excitatory, inhibitory and total conductance in a window of 50 ms before the spike 
in a cortical neuron subject to fluctuating conductance injection. The two states, equal conductances (left) 
and inhibition-dominated (right), were recreated similar to the model of Fig.|4j Conductance STAs showed 
qualitatively similar patterns. B. Vector representation showing the variation of synaptic conductances 
preceding each spike (as in Fig. |4p). 
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Figure 6: Geometrical interpretation of the average conductance patterns preceding spikes and test 
in dynamic-clamp. A. Red ellipses: isoprobable conductance configurations. Green area: con- 
ductance configurations for which the total synaptic current is positive at spike threshold. The 
vector representations of Figs. |^-[5] are schematized here, and compared to the lines defined by 
{ge{Ee - Vt) +gi{Ei - Vt) + Gi{Ei - Vt) = 0} (solid gray line) and {ge+gi = geO+gio} (dashed gray line). 
The angle between the two lines and the aspect ratio of the ellipse determine whether spikes are preceded, 
on average, by total conductance increase (left) or decrease (right) (see text for further explanations). B. 
Spike-triggered average conductances obtained in dynamic-clamp, illustrating that for the same average 
conductances, the variances determine whether spikes are preceded by total conductance increase (left) or 
decrease (right). C. Geometrical prediction tested in dynamic-clamp (left): grouped data showing total 
conductance change preceding spikes as a function of the ratio ag/o,. The dashed line {Oe/Oi=0.6) visual- 
izes the predicted value separating total conductance increase cases from total conductance decrease cases. 
In addition (right), dynamic-clamp data indicates that the amplitude of change of each of the conductances 
before a spike is linearly correlated with the standard deviation parameter used for this conductance. 
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Figure 7: Spike-triggered conductance extraction from intracellular recordings of the V,,, and test in 
dynamic-clamp. A. Left: spike-triggered average (STA) of the V,„ in an integrate-and-fire extension of 
the point-conductance model (top trace). The numerically obtained conductance STAs (orange, green) are 
compared to the conductance STAs extracted from the V,„ (black) (bottom trace). Right: test of the STA 
method using dynamic-clamp. The STA of the Vm is obtained following injection of fluctuating conduc- 
tances (top trace). The measured conductance STAs (orange, green) are compared to the conductance STAs 
extracted from the (black) (bottom trace). B. Grouped data comparing conductance STAs extracted us- 
ing the method with the conductance STAs measured following dynamic-clamp injection: amplitude of 
conductance change preceding the spike (top graphs) and time constant of this change (bottom graphs), 
for both excitation and inhibition. C. Top: correlation between errors for excitation and inhibition on the 
absolute value of conductance variation. Bottom: total conductance change preceding spikes; comparison 
between extracted and measured STAs. Dashed lines: Y = X. 
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Figure 8: Analysis of the average conductance patterns preceding spikes, same analysis as Fig. [6^-C, but 
using conductance STAs extracted from the V^, instead of the measured ones. A. Example conductance 
STAs extracted from the V,„ STAs of the same cell and the same conductance injections as Fig. |6p. B. 
Left: test of the geometrical prediction (dashed line) using conductance STAs extracted from the V^, and 
showing total conductance change preceding spikes as a function of the ratio Oe/o, (as in Fig. [6p, left). 
Right: correlation between the amplitude of change of each conductance preceding a spike, as extracted 
from the V,„, and the standard deviation parameter for this conductance (compare to Fig.[6p, right). 
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Figure 9: Spike-triggered conductance analysis in vivo. A. STA conductance analysis from intracellular 
recordings in awake and sleeping cats. Two example cells are shown during wakefulness, and for each, 
the Ym STA (top) and the extracted conductance STAs (bottom) are shown. In the first cell (left), the 
total conductance increases before the spike. In the second example cell (right), the total conductance 
decreases before the spike (black traces are exponential fits to the extracted STAs). B. Total conductance 
change preceding spikes as a function of the ratio Og/O/. Given the cell-to-cell variability of observed spike 
thresholds, each cell has a different predicted ratio separating total conductance increase cases from total 
conductance decrease cases. The two dashed lines (Oe/o,=0.48 and Og/o,-!. 07) visualize the two extreme 
predicted ratios. Cells in white are the ones not conforming to the prediction. C. Dependency of the ratio 
Oe/O; estimated by the VmD method on the value of the membrane capacitance C. Two sets of realistic 
Vm distribution parameters were used as input for the estimation, one leading to Og/o, > 0.6 (left), another 
leading to Og/O/ < 0.6 (right). For each set, the total input resistance was varied from 10 MQ. (bottom 
curves) to 50 MQ. (top curves), in steps of 10 MQ.. Panel A modified from Rudolph et al., 2007. 



